翻訳と辞書
Words near each other
・ "O" Is for Outlaw
・ "O"-Jung.Ban.Hap.
・ "Ode-to-Napoleon" hexachord
・ "Oh Yeah!" Live
・ "Our Contemporary" regional art exhibition (Leningrad, 1975)
・ "P" Is for Peril
・ "Pimpernel" Smith
・ "Polish death camp" controversy
・ "Pro knigi" ("About books")
・ "Prosopa" Greek Television Awards
・ "Pussy Cats" Starring the Walkmen
・ "Q" Is for Quarry
・ "R" Is for Ricochet
・ "R" The King (2016 film)
・ "Rags" Ragland
・ ! (album)
・ ! (disambiguation)
・ !!
・ !!!
・ !!! (album)
・ !!Destroy-Oh-Boy!!
・ !Action Pact!
・ !Arriba! La Pachanga
・ !Hero
・ !Hero (album)
・ !Kung language
・ !Oka Tokat
・ !PAUS3
・ !T.O.O.H.!
・ !Women Art Revolution


Dictionary Lists
翻訳と辞書 辞書検索 [ 開発暫定版 ]
スポンサード リンク

Rayleigh-Taylor instability : ウィキペディア英語版
Rayleigh–Taylor instability

The Rayleigh–Taylor instability, or RT instability (after Lord Rayleigh and G. I. Taylor), is an instability of an interface between two fluids of different densities which occurs when the lighter fluid is pushing the heavier fluid.〔
Drazin (2002) pp. 50–51.〕 Examples include the behavior of water suspended above oil in the gravity of Earth,〔 mushroom clouds like those from volcanic eruptions and atmospheric nuclear explosions,〔http://gizmodo.com/why-nuclear-bombs-create-mushroom-clouds-1468107869〕 supernova explosions in which expanding core gas is accelerated into denser shell gas,〔. See page 274.〕
and instabilities in plasma fusion reactors.
Water suspended atop oil is an everyday example of Rayleigh–Taylor instability, and it may be modeled by two completely plane-parallel layers of immiscible fluid, the more dense on top of the less dense one and both subject to the Earth's gravity. The equilibrium here is unstable to any perturbations or disturbances of the interface: if a parcel of heavier fluid is displaced downward with an equal volume of lighter fluid displaced upwards, the potential energy of the configuration is lower than the initial state. Thus the disturbance will grow and lead to a further release of potential energy, as the more dense material moves down under the (effective) gravitational field, and the less dense material is further displaced upwards. This was the set-up as studied by Lord Rayleigh.〔 The important insight by G. I. Taylor was his realisation that this situation is equivalent to the situation when the fluids are accelerated, with the less dense fluid accelerating into the more dense fluid.〔 This occurs deep underwater on the surface of an expanding bubble and in a nuclear explosion.
As the RT instability develops, the initial perturbations progress from a linear growth phase into a non-linear or "exponential" growth phase, eventually developing "plumes" flowing upwards (in the gravitational buoyancy sense) and "spikes" falling downwards. In general, the density disparity between the fluids determines the structure of the subsequent non-linear RT instability flows (assuming other variables such as surface tension and viscosity are negligible here). The difference in the fluid densities divided by their sum is defined as the Atwood number, A. For A close to 0, RT instability flows take the form of symmetric "fingers" of fluid; for A close to 1, the much lighter fluid "below" the heavier fluid takes the form of larger bubble-like plumes.〔
This process is evident not only in many terrestrial examples, from salt domes to weather inversions, but also in astrophysics and electrohydrodynamics. RT instability structure is also evident in the Crab Nebula, in which the expanding pulsar wind nebula powered by the Crab pulsar is sweeping up ejected material from the supernova explosion 1000 years ago. The RT instability has also recently been discovered in the Sun's outer atmosphere, or solar corona, when a relatively dense solar prominence overlies a less dense plasma bubble. This latter case is an clear example of the magnetically modulated RT instability.〔. See Chap. X.〕
Note that the RT instability is not to be confused with the Plateau–Rayleigh instability (also known as Rayleigh instability) of a liquid jet. This instability, sometimes called the hosepipe (or firehose) instability, occurs due to surface tension, which acts to break a cylindrical jet into a stream of droplets having the same volume but lower surface area.
Many people have witnessed the RT instability by looking at a lava lamp, although some might claim this is more accurately described as an example of Rayleigh–Bénard convection due to the active heating of the fluid layer at the bottom of the lamp.
==Linear stability analysis==

The inviscid two-dimensional Rayleigh–Taylor (RT) instability provides an excellent springboard into the mathematical study of stability because of the simple nature of the base state.〔Drazin (2002) pp. 48–52.〕 This is the equilibrium state that exists before any perturbation is added to the system, and is described by the mean velocity field U(x,z)=W(x,z)=0,\, where the gravitational field is \textbf=-g\hat\quad \gamma=} \quad\text\quad \mathcal=\frac}}}},\,
where \gamma\, is the temporal growth rate, \alpha\, is the spatial wavenumber and \mathcal\, is the Atwood number.
The perturbation introduced to the system is described by a velocity field of infinitesimally small amplitude, (u'(x,z,t),w'(x,z,t)).\, Because the fluid is assumed incompressible, this velocity field has the streamfunction representation
:\textbf'=(u'(x,z,t),w'(x,z,t))=(\psi_z,-\psi_x),\,
where the subscripts indicate partial derivatives. Moreover, in an initially stationary incompressible fluid, there is no vorticity, and the fluid stays irrotational, hence \nabla\times\textbf'=0\,. In the streamfunction representation, \nabla^2\psi=0.\, Next, because of the translational invariance of the system in the ''x''-direction, it is possible to make the ansatz
:\psi\left(x,z,t\right)=e^\Psi\left(z\right),\,
where \alpha\, is a spatial wavenumber. Thus, the problem reduces to solving the equation
:\left(D^2-\alpha^2\right)\Psi_j=0,\,\,\,\ D=\frac,\,\,\,\ j=L,G.\,
The domain of the problem is the following: the fluid with label 'L' lives in the region -\infty, while the fluid with the label 'G' lives in the upper half-plane 0\leq z<\infty\,. To specify the solution fully, it is necessary to fix conditions at the boundaries and interface. This determines the wave speed ''c'', which in turn determines the stability properties of the system.
The first of these conditions is provided by details at the boundary. The perturbation velocities w'_i\, should satisfy a no-flux condition, so that fluid does not leak out at the boundaries z=\pm\infty.\, Thus, w_L'=0\, on z=-\infty\,, and w_G'=0\, on z=\infty\,. In terms of the streamfunction, this is
:\Psi_L\left(-\infty\right)=0,\qquad \Psi_G\left(\infty\right)=0.\,
The other three conditions are provided by details at the interface z=\eta\left(x,t\right)\,.
''Continuity of vertical velocity:'' At z=\eta, the vertical velocities match, w'_L=w'_G\,. Using the streamfunction representation, this gives
:\Psi_L\left(\eta\right)=\Psi_G\left(\eta\right).\,
Expanding about z=0\, gives
:\Psi_L\left(0\right)=\Psi_G\left(0\right)+\text,\,
where H.O.T. means 'higher-order terms'. This equation is the required interfacial condition.
''The free-surface condition:'' At the free surface z=\eta\left(x,t\right)\,, the kinematic condition holds:
:\frac+u'\frac=w'\left(\eta\right).\,
Linearizing, this is simply
:\frac=w'\left(0\right),\,
where the velocity w'\left(\eta\right)\, is linearized on to the surface z=0\,. Using the normal-mode and streamfunction representations, this condition is c \eta=\Psi\,, the second interfacial condition.
''Pressure relation across the interface:'' For the case with surface tension, the pressure difference over the interface at z=\eta is given by the Young–Laplace equation:
:p_G\left(z=\eta\right)-p_L\left(z=\eta\right)=\sigma\kappa,\,
where ''σ'' is the surface tension and ''κ'' is the curvature of the interface, which in a linear approximation is
:\kappa=\nabla^2\eta=\eta_.\,
Thus,
:p_G\left(z=\eta\right)-p_L\left(z=\eta\right)=\sigma\eta_.\,
However, this condition refers to the total pressure (base+perturbed), thus
:\left()-\left()=\sigma\eta_.\,

(As usual, The perturbed quantities can be linearized onto the surface ''z=0''.) Using hydrostatic balance, in the form
:P_L=-\rho_L g z+p_0,\qquad P_G=-\rho_G gz +p_0,\,
this becomes
:p'_G-p'_L=g\eta\left(\rho_G-\rho_L\right)+\sigma\eta_,\qquad\textz=0.\,
The perturbed pressures are evaluated in terms of streamfunctions, using the horizontal momentum equation of the linearised Euler equations for the perturbations,
: \frac = - \frac\frac\, with i=L,G,\,
to yield
:p_i'=\rho_i c D\Psi_i,\qquad i=L,G.\,
Putting this last equation and the jump condition on p'_G-p'_L together,
:c\left(\rho_G D\Psi_G-\rho_L D\Psi_L\right)=g\eta\left(\rho_G-\rho_L\right)+\sigma\eta_.\,
Substituting the second interfacial condition c\eta=\Psi\, and using the normal-mode representation, this relation becomes
:c^2\left(\rho_G D\Psi_G-\rho_L D\Psi_L\right)=g\Psi\left(\rho_G-\rho_L\right)-\sigma\alpha^2\Psi,\,
where there is no need to label \Psi\, (only its derivatives) because \Psi_L=\Psi_G\,
at z=0.\,
; Solution
Now that the model of stratified flow has been set up, the solution is at hand. The streamfunction equation \left(D^2-\alpha^2\right)\Psi_i=0,\, with the boundary conditions \Psi\left(\pm\infty\right)\, has the solution
:\Psi_L=A_L e^,\qquad \Psi_G = A_G e^.\,
The first interfacial condition states that \Psi_L=\Psi_G\, at z=0\,, which forces A_L=A_G=A.\, The third interfacial condition states that
:c^2\left(\rho_G D\Psi_G-\rho_L D\Psi_L\right)=g\Psi\left(\rho_G-\rho_L\right)-\sigma\alpha^2\Psi.\,
Plugging the solution into this equation gives the relation
:Ac^2\alpha\left(-\rho_G-\rho_L\right)=Ag\left(\rho_G-\rho_L\right)-\sigma\alpha^2A.\,
The ''A'' cancels from both sides and we are left with
:c^2=\frac\frac+\frac.\,
To understand the implications of this result in full, it is helpful to consider the case of zero surface tension. Then,
:c^2=\frac\frac,\qquad \sigma=0,\,
and clearly
* If \rho_G<\rho_L\,, c^2>0\, and ''c'' is real. This happens when the
lighter fluid sits on top;
* If \rho_G>\rho_L\,, c^2<0\, and ''c'' is purely imaginary. This happens
when the heavier fluid sits on top.
Now, when the heavier fluid sits on top, c^2<0\,, and
:c=\pm i \sqrt},\qquad \mathcal=\frac,\,
where \mathcal\, is the Atwood number. By taking the positive solution, we see that the solution has the form
:\Psi\left(x,z,t\right)=Ae^\exp\left()=A\exp\left(\alpha\sqrt}t\right)\exp\left(i\alpha
x-\alpha|z|\right)\,
and this is associated to the interface position ''η'' by: c\eta=\Psi.\, Now define B=A/c.\,
The time evolution of the free interface elevation z = \eta(x,t),\, initially at \eta(x,0)=\Re\left\,\, is given by:
:\eta=\Re\left\\,t\right)\exp\left(i\alpha x\right)\right\}\,
which grows exponentially in time. Here ''B'' is the amplitude of the initial perturbation, and \Re\left\\, denotes the real part of the complex valued expression between brackets.
In general, the condition for linear instability is that the imaginary part of the "wave speed" ''c'' be positive. Finally, restoring the surface tension makes ''c''2 less negative and is therefore stabilizing. Indeed, there is a range of short waves for which the surface tension stabilizes the system and prevents the instability forming.

抄文引用元・出典: フリー百科事典『 ウィキペディア(Wikipedia)
ウィキペディアで「Rayleigh–Taylor instability」の詳細全文を読む



スポンサード リンク
翻訳と辞書 : 翻訳のためのインターネットリソース

Copyright(C) kotoba.ne.jp 1997-2016. All Rights Reserved.